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Abstract 

We show that electron-ion structure factors in fluid metallic systems can be 
well understood from an application of linear response in the electron sys- 
tem, combined with hard-sphere like correlation for the ionic component. In 
particular, we predict that electron-ion structure factors fall into two general 
classes, one for high (Z > 3) and one for low {Z < 2) valence metals, and make 
suggestions for experiments to test these ideas. In addition, we show how the 
general success of electronic linear response for most metallic systems stems 
in part from an intrinsic interference between atomic and electronic length 
scales which weakens the nonlinear response. The main exception to this is 
metallic hydrogen. 



I. INTRODUCTION 



At a near fundamental level, liquid metals are complex binary fluids consisting of ions 
in a sea of conduction electrons, their physical properties linked to the three corresponding 
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correlation functions: Sii{k), Sei{k), and Seeik) |T[. We present a simple and evidently 
accurate analytic scheme to calculate electron-ion structure factors {Sei{k)) by combining a 
hard-sphere approximation for the ionic structure with a simple linear response theory for 
the electrons. These structure factors are now in principle accessible experimentally through 
recent advances in both neutron and x-ray scattering techniques. Another route to effective 
electron-ion interactions therefore opens, but now through the fluid state. 

We also address the evident success of the linear approximation by studying a related 
problem, the density, p*"'^(k), of an initially uniform electron gas induced by an embedded 
pseudo-potential, v^^(k). By comparing linear and second order response to full (Kohn- 
Sham 0) non-linear response, we show that even though f^'^(k) is not necessarily a small 
perturbation, the consequent response series converges term by term. The non-linear terms 
are significantly reduced by an interference between atomic and electronic length scales for 
most metals, the main exception to this being hydrogen. 

II. DISCUSSION 
A. Electron-ion correlation functions 

The electron-ion-structure factor can be written as |]T[: 

Seiik) = < /)«(k)pi'^(k) >= ^Snik), (1) 

where pf\k) is the Fourier transform of the one-particle density operator of component j, 
Sii{k) is the ion-ion structure factor and n{k) is identified as the pseudo-electron density, or 
pseudo-atom (of valence Z). Thus electron- ion correlations can be described by convolving 
the pseudo-atom with the ionic correlations. The ionic correlations are themselves well 
described by a Percus-Yevick hard-sphere structure factor 0, while the pseudo-atom is 
described by a standard linear response formulation: 

nik)=xiik)v^'{k) (2) 
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where Xi(^) is the well known linear response function; to approximate Xi(^) "we use a 
Local Density Approximation (LDA) local field factor ||^. The electron-ion interaction 
is modeled by a simple local one-parameter empty-core pseudopotential i.e.: v^'^{k) = 
— {4:7!-e^/k'^)cos{kRc), where Rc is the core radius and the pseudo-potential goes through 
zero at ko = 7t/2Rc. (We note that at this linear level, the effects of ionic averaging on 
the pseudo-atom are ignored Using the approximations above in (|l]), we compare our 
approach in Fig. 1 to the full ab-initio Car-Parrinello calculations of de Wijs a/ [§. 
The correspondence is striking, especially when we note that the parameters rj and Rc are 
a priori set by other physical properties (no fitting is necessary). 

Besides a semi-quantitative description of electron-ion structure factors, this linear re- 
sponse theory now provides an important qualitative insight into the form of the electron-ion 
structure factors 0. The pseudo-atom density, n{k), is typically largest for smaller k and 
thereafter rapidly declines for larger k, while the near classical ion-ion structure factor, 
Sii{k), follows an inverse behavior; it is small for small k. Together with the product 
form (|I|) this implies that the shape of the electron-ion structure factor, Sei{k), is deter- 
mined primarily by the the position of the zero-crossing, k^, of n{k) with respect to the first 
maximum, kp, of Sjj{k). If kQ < kp, then Sjj{k) selects (or filters) the negative part of n{k) 
and Sei{k) takes a form similar to that of Mg (Fig. 1 (a)). Conversely, if ko > kp, then the 
ion-ion structure factor selects (or filters) the positive part of n{k), and again, Sei{k) takes a 
form similar to that of Bi (Fig. 1 ([l| (c)). Since Xi(^) is positive definite, the zero-crossing 
in linear response occurs at kQ. The large slope of n{k) near the zero-crossing then implies 
that non-linear corrections must have a small effect on the location of the zero-crossing, and 
this, together with the expected accuracy of linear response, implies that fco ~ ^o- For most 
metals, k^ is just a little less than 2kF, and the latter's ratio to kp is well known: for low 



valence {Z < 2), 2kp < kp] for high valence {Z > 3): 2kp > kp |T0[. This accounts in a 
straightforward way for the two separate forms found by deWijs et a/ ^j: For Mg, ko < kp 
[Z = 2), which belongs to the low valence class of electron-ion structure factors. For 
Bi, ko > kp {Z = 5) and we may refer to this as the high valence class of electron-ion 



structure factors |TT[]. Generally ions of valence Z < 2 belong to the low valence class while 
ions with valence Z > 3 belong to the high valence class. Ions with valence Z = 3 typically 
belong to the high valence class also, although they may be characterized by a crossover 
form . The analytical approach above can easily be extended by using the modern theory 
of classical liquids to obtain improved ion- ion structure factors |T2|, but to include second 
order contributions to the pseudo-atom n{k) necessitates not only second order electron re- 
sponse, but also contributions from ion-ion triplet structure. The latter can be carried out 
with concepts from the theory of classical liquids p. 

B. Proposed Experiments 

The principal features of electron-ion structure factors can be measured by exploiting 
the differences between the x-ray scattering structure factor, Sfj{k), determined with a free- 
atom form factor, fA^k), and the structure factor, Sf^j{k), determined by neutron scattering 



||T3| , p!^ . As emphasized by Chihara |]14| the x-ray structure factor for liquid metals will equal 
Sf^j{k) only when determined with an ionic form factor augmented by a pseudo-atom form 
factor; i.e. fi{k) + n{k), so that: 



SMk) \fi{k)+n{k)\' 



(3) 
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This effect is clearly expected to be largest for metals with larger ratios of valence to core 
electrons. Thus we predict a small effect for metals with smaller valence to core ratio such 
as Na or K, a 2% difference at the 1st peak of the structure for Li (ratio= 1 : 2) or Al 
(ratio= 3 : 10), but by far the largest effect for Be (ratio = 1:1) where the difference at 
the principal peak of the structure factor could be as high as 7%, well within experimental 
range. Another interesting candidate would be metallic Si (ratio = 4 : 10) since covalent 
effects still make themselves felt in the liquid state suggesting that experiments could reveal 
effects beyond linear response. 

To date the experimental electron-ion structure factors and related pseudo-atoms show 
considerably more structure than indicated by theoretical predictions |T^. Significant ex- 



perimental challenges are faced in the accuracy resulting from subtraction of two sets of 
data obtained by quite different means, each with important (but different) systematic cor- 
rections; however the present approach suggests that the current differences between x-ray 



and neutron scattering should be reexamined (see also Jl^)- Using a pseudo-atom instead of 
the full free atom as a form factor can assist in comparing neutron and x-ray measurements 
and help unravel various systematic corrections applied. The advent of high precision x-ray 
and neutron sources currently coming on-line suggests that these proposed effects can be 
systematically explored. 



C. Non-linear response of an atom in an electron gas 

The evident (and long-standing) success of the linear response approximation for electron 
response here demonstrated for electron- ion structure factors, calls for further investi- 
gation. The accuracy of linear response in a crystalline solid is commonly attributed to the 
fact that the structure dependent reciprocal lattice vectors are typically near the pseudo- 
potential zero-crossing, ko, with the associated inference that the net scattering is smaller 
than one would naively expect @]. For liquids or other disordered systems such arguments 
are less appropriate. To examine the strength of linear response in the absence of ionic struc- 
ture, we consider a simpler problem, namely the response of the interacting electron gas to 
a single ion, where the electron-ion interaction is modeled by a simple local one-parameter 
empty-core pseudopotential P]. We will compare two routes to the induced density, p*"'^(fc). 
The first follows from solving the Kohn-Sham equations [0 exactly (for the given pseu- 
dopotential) within the local density approximation (LDA), the second from the standard 
expansion of the response in powers of the perturbing (pseudo)potential, i.e.; 

r\k) = XiikV^k) X2{k, fci, k2)v^\k,)v''%k2) + ■■■. (4) 

ki k2 

Here the response functions, Xn{ki---), are properties of the homogeneous interacting electron 
gas. In Fig. 2 we compare the explicit second order response with an LDA local field 
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factor P,p!6| to the full non-linear LDA response. Clearly the non-linear response is well 



characterized by the 2nd order term, implying that the success of linear response is not due 
to cancellation between higher order terms of opposite sign but instead that each successive 
term is individually small compared to the previous term in the expansion; the response 
series converges very rapidly, term by term. 

The non-linear response is largest for atomic parameter Rc = (hydrogen), and decreases 
with a larger atomic-parameter, Rc, as might be physically anticipated. However as Rc 
increases from zero, a noticeable secondary minimum occurs when the inverse atomic length, 
ko, is equal to 2/ci?. For the cases plotted in Fig. 2, the maximum in second order response 
at kQ/2kp = 1 (or Rc/r^ = 0.41) is reduced by an entire order of magnitude when compared 
with the maximum in second order response calculated for hydrogen {Rc/vs = 0), and is 
typically equal to the value at SRc/vg. The physics behind this minimum is attributed to 
the following; the second order response function, X2{k,ki,k2), peaks when the summed 
arguments in (||) are close to 2kp 0. If the pseudo-potential zero-crossing, ko, is near the 
response peaks at 2k p, a maximal cancellation or maximal destructive interference of the 
atomic and electronic length scales occurs, leading to a minimum in second order response. 
The ratio of the atomic and electronic length scales is set primarily by the volume energy 
terms in the total ground state energy, and is almost independent of structure [Q; kQ/2kp lies 
between 0.75 and 1 for most metals, and is therefore very close to the secondary minimum 
in the non-linear response. 

As noted, the effect we discuss originates from an interference between intrinsic atomic 
and electronic length scales, but it also complements the argument given for crystalline 
solids alone, which stems from the confluence of an atomic and a structural length scale. 
The clear exception to these interference effects is again the singular case of a point-charge 
{v^^{k) ~ 47re^/fc^), i.e. the case of hydrogen, which has no well-defined core-length scale, 
ko, no oscillations in the potential and thus no interference effect in the higher order terms. 
In contrast to other systems, non-linear response terms are large term by term. In fact, 
the response series may not even formally converge and care must be taken when applying 
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concepts derived from linear-response theory to hydrogen (it is not a simple material) . 

Finally we note that the second order response contribution is of the same order as the 
difference between first order response with or without local field corrections. In addition, 
the combined effects of exchange and correlation partially cancel between first and second 
order, implying that neglect of higher order response results in an over- estimation of the role 
of exchange and correlation, which, in turn, has important implications for the widespread 
application of linear response theory in the derivation of effective ion-ion potentials in (sim- 
ple) metals. 

III. CONCLUSIONS 

The electron-ion structure factors of liquid simple-metals are well described by a sim- 
ple linear response theory augmented by linear response for the electrons. This approach 
suggests two main classes of electron-ion correlation functions, one for high and one for 
low valence metals. Experimental advances in x-ray and neutron-scattering may be able to 
provide measurements of these electron-ion correlation functions, with liquid Be being the 
most promising candidate. A route to information on fundamental electron- ion interactions 
therefore becomes available through the fluid state. Finally, the well documented success 
of the linear response approximation for electrons stems in part from an interference effect 
between atomic and electronic length scales. 

This work was supported by the NSF through the Cornell Center for Materials Research 
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FIGURE CAPTIONS 



FIGURE 1 

The electron-ion structure factors Sei{k) and related electron- ion correlation functions 
Qeiir) for Mg and Bi: Car-Parrinello results of de Wijs et al (solid line) vs. the simple 
linear-response approach augmented by a hard-sphere approximation (dashed line). Panel 
(a) shows Seiik) and panel (b) shows Qeiij) for liquid Mg. Panel (c) shows Seiik) and 
panel (d) shows gei{f) for liquid Bi. For Mg the parameters (taken from the literature) are: 
rs=2.66ao, Rc=1.31aQ and for Bi the parameters (taken from the literature) are: rs=2.25ao 
and Rc=1.15aQ. Both have a hard-sphere parameter, ri=0A6. (note that for the 5'e/(^) the 
region inside the core radius is not physically significant.) 



FIGURE 2 

A comparison of full non-linear LDA response [p(k) — p^^\k)] (solid line) to second order 
LDA response (dashed line) for an empty core pseudo-potential with Rc=l.5ao embedded 
in an electron gas with density parameter, rs=3aQ. For the scale, compare this to the full 
response with the limit p{k — > 0)=1. The higher order response is of the order of a few % of 
the full response and in turn, the second order response captures almost all the non-linear 
response. (The small difference at A; — is a numerical artifact stemming from the use of 
a large but finite real-space cut-off radius in the Kohn-Sham procedure.) In the insert is 
plotted the maximum of the 2nd order response vs. Rc/ts for rs=2ao(dotted), rs=3ao(solid) 
and rs=5ao (dashed). Note especially the minimum at Rc/rs=OAl which corresponds to 
ko=2kF- It is reduced by an order of magnitude from the value at Rc=0 (hydrogen) and is 
traced to an interference between atomic and electronic length-scales. 
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